Attorney's Docket No.: 07844-647001 / P600 



APPLICATION 



FOR 



UNITED STATES LETTERS PATENT 



TITLE: CONSTRUCTING COLOR PROFILE FOR 

PARAMETERIZED IMAGE COLOR SPACE 



APPLICANT: MICHAEL BOURGOIN AND LARS U. BORG 



CERTIFICATE OF MAILING BY EXPRESS MAIL 
Express Mail Label No. EV3993121 15US 



Date of Deposit 



April 7. 2004 



Patent Application 

Att Dkt. 07844-647001 / P600 



CONSTRUCTING COLOR PROFILE FOR PARAMETERIZED IMAGE 

COLOR SPACE 

BACKGROUND OF THE INVENTION 
[0001] The present application describes systems and techniques relating to color 
management, for example, supporting color management of an image with a parameterized 
image color space. 

[0002] Devices that are capable of representing image data (monitors, printers, scanners, 
cameras, etc.) often have significant variations in color gamut, the range of colors producible by 
the device. To accommodate this range in device gamuts, a variety of color spaces, and color 
management systems and techniques have been developed. Color management enables 
different color space values to be mapped from one device's gamut to another's using color 
profiles, which define how color information is transformed into or out of a standard reference 
space called a profile connection space. 

[0003] Using color management, a human's perception of color from one device 
representation to another can be kept close to constant, despite the variation in color capabilities 
of different devices. To assist in color management, the Commission Internationale de 
L'Eclairage (CIE) has defined various color spaces that are device-independent and encompass 
the full gamut of human vision, and can thus be used as profile connection spaces. Typical 
profile connection spaces in color management systems include CIEXYZ and CIELAB. LAB 
space is a color space having a luminance channel, L, and opponent color channels, A 
(green<->red) and B (blue«*yellow). 

[0004] A color profile defines how to transform color information from one color space to 
another, such as from a device-dependent color space into a profile connection space, or the 
reverse. Many color profiles also conform to a defined color profile architecture, which 
provides flexibility in their use. For example, the International Color Consortium (ICC) 
provides a defined color profile architecture commonly used in many color management 
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systems. ICC profiles have been developed for many different color spaces, which typically 
have standardized component ranges to facilitate their use. 

SUMMARY OF THE INVENTION 
[0005] In general, the invention features systems and techniques relating to supporting color 
management of an image with a parameterized image color space. According to an aspect, a 
color profile that conforms to a defined color profile architecture and that defines a multistage 
transform capable of translating a first color space to a second color space is generated. An 
image can be processed using the color profile, and the image can be output. The image 
includes a parameterized encoding of an image color space with image parameters defining a 
range and an offset of an image component of the image, and a white point of the image color 
space. One of the first and second color spaces can be the image color space while the other of 
the first and second color spaces can be a profile connection space. Moreover, generation of the 
color profile involves affecting two stages of the multistage transform based on the image 
parameters. 

[0006] Generating the color profile can involve increasing image processing precision by 
fitting output to input data scopes between the two stages of the multistage transform based on 
the parameterized encoding of the image. Generating the color profile can involve affecting 
three stages, which include the two stages, of the multistage transform such that the color 
profile effects chromatic adaptation according to the white point and transcodes the image 
component according to the range and the offset, where the fitting includes fitting output to 
input data scopes among the three stages. The image parameters of the parameterized encoding 
can define ranges, offsets, and bit depths of image components of the image, and the color 
profile can be a bit-depth independent color profile. 

[0007] Affecting the two stages can involve determining a first processing stage of a 
transform-defining element of the defined color profile architecture, wherein the first 
processing stage accounts for the range and the offset; determining a second processing stage of 
the transform-defining element, wherein the second processing stage defines at least a portion 
of a conversion of the image color space to a chromatic adaptation color space; and determining 
a third processing stage of the transform-defining element, wherein the third processing stage 
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defines a chromatic adaptation in the chromatic adaptation color space according to the white 
point. Affecting the two stages can involve determining a first processing stage of a transform- 
defining element of the defined color profile architecture, wherein the first processing stage 
accounts for the range; determining a second processing stage of the transform-defining 
element, wherein the second processing stage accounts for the offset and defines at least a 
portion of a conversion of the image color space to a chromatic adaptation color space; and 
determining a third processing stage of the transform-defining element, wherein the third 
processing stage defines a chromatic adaptation in the chromatic adaptation color space 
according to the white point. 

[0008] According to another aspect, a system includes a device, and a data processing 
machine including an input-output interface, an operating system, and a color management 
software component that generates a color profile for an image. The image includes a 
parameterized encoding of an image color space with image parameters defining ranges and 
offsets of image components of the image. The color management software component 
generates the color profile by representing a transformation including a matrix followed by 
curves in a pipeline of a defined color profile architecture while increasing processing precision 
governed by the color profile based on the parameterized encoding. The image parameters can 
further define a white point of the image color space, and the transformation can include the 
matrix followed by the curves followed by an additional matrix. The color profile can effect 
chromatic adaptation according to the white point and can transcode the image components 
according to the ranges and the offsets, in at least three stages of the pipeline. The color profile 
can transcode the image components in a multidimensional interpolation table stage and a one 
dimensional transforms stage of the pipeline, and the color profile can effect chromatic 
adaptation in a matrix stage of the pipeline. 

[0009] According to another aspect, a machine-implemented technique includes improving 
color accuracy of conversion of an image color space of an image by affecting two or more 
processing stage definitions for a defined image processing pipeline, based on image 
parameters, such that the defined image processing pipeline transcodes an image component 
according to a range and an offset. The affecting includes fitting output to input data scopes 
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between two of the processing stage definitions, and the image includes a parameterized 
encoding of the image color space with the image parameters defining the range and the offset 
of the image component of the image. The image parameters can further define a white point of 
the image color space, and the affecting can involve affecting three or more processing stage 
definitions for the defined image processing pipeline, based on the image parameters, such that 
the defined image processing pipeline transcodes the image component according to the range 
and the offset, and effects chromatic adaptation according to the white point. 
[0010] The invention can be implemented to realize one or more of the following 
advantages. The systems and techniques described here involve construction of color profiles 
that accommodate nonstandard data encodings for a color space. The color profiles constructed 
can be independent of the color space's data bit depth. Moreover, the color profiles constructed 
can handle transcoding (decoding and/or encoding) of the nonstandard data encodings and 
provide improved color accuracy. In one embodiment, an ICC profile is constructed that 
transcodes, color converts and performs chromatic adaptation on a parameterized encoding of 
an image color space, while also maximizing precision and preventing clipping of valid data. 
This can prevent color banding and color shifts that might otherwise occur. By distributing 
parts of the image processing through the stages of a defined processing pipeline based on the 
encoding of an image color space, quantization error and color error can be minimized. 
[0011] The details of one or more embodiments of the invention are set forth in the 
accompanying drawings and the description below. Other features and advantages of the 
invention will become apparent from the description, the drawings, and the claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0012] FIG. 1 illustrates a system with parameterized image color space color management. 
[0013] FIG. 2 illustrates a portion of a color-managed workflow in the system of FIG. 1. 
[0014] FIG. 3 illustrates example color spaces. 

[0015] FIG. 4 illustrates creation of a transform-defining element based on image 
parameters. 

[0016] FIG. 5 illustrates fitting output to input data scopes between two stages of a defined 
image processing pipeline. 
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[0017] FIG. 6 illustrates creation of a transform-defining element according to an 
implementation. 

[0018] FIG. 7 illustrates a transform-defining element according to an implementation. 

[0019] FIG. 8 shows equations according to an implementation. 

[0020] FIG. 9 illustrates creation of a transform-defining element according to another 

implementation. 

[0021] Like reference symbols in the various drawings indicate like elements. 

DETAILED DESCRIPTION 
[0022] FIG. 1 illustrates a system with parameterized image color space color management. 
A data processing machine 100 includes an input-output interface 110 that provides 
connectivity to one or more devices 150. The device(s) 150 include one or more image 
resources 160, which have nonstandard color space encodings. The device(s) 150 can be of 
various types, including storage devices (e.g., hard disks), input-output devices for the machine 
100 (e.g., a printer, monitor, scanner, etc.), or consumer electronic devices (e.g., a digital 
camera, mobile phone, portable electronic poster, digital picture frame, etc.). The device(s) 150 
can be local or remote, and can be connected to the machine 100 via a network, wireless link, 
cable, or system bus. For example, the machine and device 100, 150 can be a personal 
computer; or the device 150 can be integrated with the machine 100 into a composite device, 
such as a printer/scanner that includes a processor and performs color management operations, 
or a portable digital assistant (PDA). 

[0023] The data processing machine 100 can include an operating system 120, which is 
software used to control operation of the machine 100 and to direct the processing of programs. 
The data processing machine 100 includes a color management software component 130 that 
generates a color profile for an image having a parameterized encoding of an image color space. 
For example, an image 160 can be a parameterized CIELAB image (e.g., a Joint Photographic 
Experts Group (JPEG) image with parameterized CIELAB color space, which may have nearly 
arbitrary component ranges, offsets and bit depths, as well as a nearly arbitrary white point). 
The color management component 130 can generate a color profile capable of translating the 
parameterized CIELAB image into a profile connection space. 
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[0024] The color management component 130 can generate a color profile with decoding 
and/or encoding capability: the decoding capability being decoding a parameterized image color 
space while translating the color space into a profile connection space, and the encoding 
capability being encoding a parameterized image color space while translating a profile 
connection space into the color space. The color management component 130 can be a 
software application or plug-in for a software application, or a color matching module or color 
engine. The color management component 130 can be integrated into the operating system 120, 
or alternatively, no operating system 120 need be present. Various other implementations are 
also possible, for example, the color management component 130 can be a portion of a 
document (e.g., a container in a Portable Document Format (PDF) file), or the color 
management component 130 can be a hardware driver, such as a graphics hardware driver. 
[0025] FIG. 2 illustrates a portion of a color-managed workflow in the system of FIG. 1. A 
color profile 200 has been created by the color management component 130 and conforms to a 
defined color profile architecture 210. For example, the color profile 200 may be an ICC 
profile. The color profile 200 defines a multistage transform capable of translating from a first 
color space to a second color space, such as from a color space 220 to a profile connection 
space 230. The color profile 200 may also define an additional transform capable of translating 
from the second color space to a third color space. For example, the color profile 200 may 
translate a color space 220 into a profile connection space 230 and then into a color space 240, 
which may be a working color space of a software application. 

[0026] The color space 220 can be a parameterized encoding of an image color space for an 
image with image parameters defining a range and an offset of an image component, plus a 
white point of the image color space. For example, the color space 220 can be a parameterized 
CIELAB image color space. Conventional CDELAB color has canonical component ranges; it 
is common for applications to assume that any CIELAB data encountered uses the canonical 
ranges. In the case of ICC 8-bit per component CIELAB data, the canonical component ranges 
are as follows: L on [0.0 .. 100.0], A on [-128.0 .. 127.0], and B on [-128.0 .. 127.0]. In 
addition, the canonical white point of CIELAB color is the CIE illuminant D50 [X=0.9642, 
Y=1.0000, Z=0.8249]. 
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[0027] The assumption of canonical ranges on the part of applications and color 
management systems results in incorrect processing of the color data if it actually represents 
CIELAB with nonstandard ranges. CEELAB color data is usually encoded as 8 or 16 bit 
unsigned integer components. The CIELAB component ranges are mapped onto [0 .. 255] in 
the 8 bits per component case, or onto [0 .. 65535] in the 16 bits per component case. This 
mapping is computed using the concept of data encoding ranges and offsets: the encoding range 
for a component is the number of units in the span of the data range. For L on [0.0 .. 100.0], the 
encoding range, RL, is 100; for A and B on [-128.0 .. 127.0], the encoding ranges, RA and RB, 
are both 255. The encoding offsets are the amounts that are to be added to the minimum data 
range value to equal zero. For L, the offset, OL, is 0. For A and B on [-128.0 .. 127.0], the 
encoding offsets are both 128. Given RL, OL, RA, OA, RB, and OB, the encoding to n bits per 
component is computed as: 

_ (Z + OlX2'-l) . ^ + 04r-l) . _ (B + OBj2'-l) 
UJ enc RL ' enc ~ RA 9 enc ~ RB 

In the case of JPEG parameterized CIELAB, the default parameterization is the encoding 
specified in ITU-T.42, which is a color fax standard. The ITU-T.42 CIELAB ranges are: L on 
[0.0 .. 100.0], A on [-85.0 .. 85.0], B on [-75.0 .. 125.0]. Thus, an assumption of canonical 
ranges for CIELAB data in the case of the default ranges for JPEG 2000 parameterized JPX 
CIELAB would produce serious color errors. 

[0028] The machine 100 can generate color profiles for images that have a parameterized 
encoding of an image color space with image parameters defining a range and an offset of an 
image component of the image, and a white point of the image color space. The machine 100 
can construct a new color profile for each new nonstandard range encoding encountered, or the 
machine 100 can save previously constructed color profiles for use with other images that have 
the same parameterized ranges, offsets and white point, but different bit depths, as the color 
profiles constructed can be bit-depth independent (e.g., when the received offsets are expressed 
in the bit-depth of the incoming data, these offsets can be translated into floating point values 
and put into the color profile being built without the bit depth dependency). These color 
profiles can then be used in processing the images. For example, a color profile 200 can be 
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embedded in an image 160 (either before or after compression, such as in the case of JPEG) and 
saved to a storage device 150; or a color profile 200 can be used in a defined image processing 
pipeline to transform the image from the image color space to a working color space, and the 
transformed image can be saved to a memory device 150. Moreover, when generating the color 
profiles, the color space specified in a header of the color profile may be three color (3CLR) so 
that a color engine will be sure to use the color profile in place of another color profile to which 
the color engine may otherwise default. 

[0029] FIG. 3 illustrates example color spaces. Objects with different spectral reflectivities 
can have the same color appearance, or tristimulus values. XYZ tristimulus values form the 
basis of the CLE 1931 chromaticity diagram illustrated in FIG. 3, where x, y, and z tristimulus 
values have been normalized to remove brightness and show chromaticity in only two 
coordinates, x and y. In this diagram, a horseshoe shaped color space 300 illustrates the full 
gamut of human vision, where the spectrum locus starts at four hundred and twenty nanometers 
(nm) at a point 310 and wraps around the horseshoe to six hundred and eighty nmata point 
320. 

[0030] A color space 350 represents a standardized color space, with canonical encoding 
ranges, offsets, bit depths and a white point 360 (e.g., the CIE illuminant D50). As can be seen 
in FIG. 3, a color space can include color values that lie outside the color space 300, and are 
thus not perceivable by a human. An image color space 370 is of the same general type as the 
color space 350, but is defined by image parameters that allow for arbitrarily defined 
component ranges, offsets, bit depths, and white point. Thus, an image can have a color space 
defined uniquely for that image, including an image specific white point 380. 
[0031] FIG. 4 illustrates creation of a transform-defining element 450 based on image 
parameters 410. The image parameters 410 are additional to the general type of color space 
(e.g., LAB) and thus create an image color space, a subspace of the general color space. The 
image parameters 410 can be part of the image 400 (e.g., enumerated parameters specified in a 
field in the header of the image 400), or retained elsewhere. In general, multiple stages of the 
transform-defining element 450 are affected, based on the image parameters 410, during 
generation of the color profile. For example, the image parameters 410 can be taken into 
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account within three processing stage definitions 460, 470, 480 of the transform-defining 
element 450, which may have only three stages or more than three stages. The transform- 
defining element 450 can be a lutAtoB tag or a lutBtoA tag in an ICC profile. Example 
implementations are described in detail below for constructing a lutAtoB tag in an ICC profile 
(e.g., using the interior processing stages of the lutAtoB tag and setting the A curves and B 
curves to identity), but the systems and techniques described herein are applicable in different 
variations, including application to other transform-defining elements and other defined color 
profile architectures. 

[0032] The transform-defining element 450 represents a multistage transform of the color 
profile. Once the image parameters 410 are taken into account, the multistage transform can 
translate from a first color space to a second color space while also effecting chromatic 
adaptation according to a specified white point, and transcoding image components according 
to specified ranges, offsets and bit depths. Moreover, the process of taking the image 
parameters 410 into account during generation of the color profile can involve fitting output to 
input data scopes between two or more stages of the multistage transform based on the 
parameterized encoding of the image, increasing image processing precision and avoiding data 
value clipping. 

[0033] FIG. 5 illustrates fitting output to input data scopes between two stages of a defined 
image processing pipeline 500. The processing pipeline 500 includes stages 1-5 (e.g., a 
pipeline defined by an ICC lutAtoBtag). An image 550, having a parameterized image color 
space, is to be translated by the pipeline 500 into an image 560 in a profile connection space 
(PCS). Decoding of the parameterized image color space can be performed in two stages 510, 
520 of the defined image processing pipeline 500, and the way in which the image parameters 
are taken into account within these two stages can affect the quality of the resulting image data 
in the image 560. 

[0034] For example, in a first color profile, an output data scope 5 12 of the stage 5 10 for 
the image 550 is larger than an input data scope 522 of the stage 520. This results in clipping of 
image data. In a second color profile, an output data scope 514 of the stage 510 for the image 
550 is smaller than an input data scope 524 of the stage 520. This results in loss of precision. 
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In contrast, in a third color profile, an output data scope 516 of the stage 510 for the image 550 
has been fit to an input data scope 526 of the stage 520. This results in increased precision and 
improved color accuracy during image processing in the pipeline 500. 

[0035] Spreading out the application of color space decoding parameters in a defined image 
processing pipeline can help in fitting output to input data scopes, as the machine representation 
of the values defining the stages (e.g., the representation of curves and matrix values in the 
color profile) has a limited precision in the expression of those values (e.g., a 16 bit precision 
limit for expressing fractions). The range of values that can be pushed through the pipeline is 
fit to the image color space to be processed using the color profile. The range of theoretical 
values that could be processed by a pipeline stage without clipping, based on the color profile 
generated, is the same as the range of actual input values to be fed into the pipeline stage, based 
on the image color space. Thus, the quantization error, expressed in terms of the output of a 
pipeline stage, varies with the range of input values possible in the given image color space. 
Fitting of output to input data scopes can be applied to two or more stages, or all stages, of a 
transform-defining element of a color profile. Example techniques for performing fitting of 
output to input data scopes are described below. 

[0036] FIG. 6 illustrates creation of a transform-defining element according to an 
implementation. A first processing stage that accounts for the range and the offset of a 
parameterized color space is determined at 600. A second processing stage that defines at least 
a portion of a conversion of the image color space to a chromatic adaptation color space is 
determined at 610. A third processing stage that defines a chromatic adaptation according to 
the white point of the parameterized color space is determined at 620. 

[0037] FIG. 7 illustrates a transform-defining element according to an implementation. The 
first processing stage determined above can be a multidimensional interpolation table 760. The 
second processing stage determined above can be a set of one dimensional transforms 770 (e.g., 
M curves in the lutAtoBtag). The third processing stage determined above can be a matrix 780. 
Thus, for example, mapping of non-canonical ranges, offsets and channel bit depths can be 
defined in a three dimensional interpolation table stage of the transform-defining element 450, 
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and non-canonical white points can be compensated for in a matrix stage of the transform- 
defining element 450. 

[0038] In the context of constructing an ICC v4 profile with a CIEXYZ PCS, the invention 
can combine the conversion from a non-canonical CIELAB encoding with the conversion from 
CIELAB to the CIEXYZ D50 PCS, as follows, resulting in color profiles that are independent 
of the data channel bit depths. However, in some color management systems, the 
representation of incoming data may be limited to a specified set of bit depths (e.g., limited to 8 
bit or 16 bit incoming data). In such cases, the system can be modified to use higher bit depths 
for data with a bit depth not in the specified set, and appropriate scaling can be applied. For 
example, in a 8 and 16 bit system, when 12 bit input is provided, a client in the system can put 
this into a 16 bit field and then apply scaling to make the data appear to be 16 bit data. 
[0039] In an ICC v4 profile with a CIEXYZ PCS, the lutAtoBtags can be constructed to 
implement a matrix/offset/ ID Luts form of the CIELAB to XYZ conversion, with mapping for 
the non-canonical ranges, offsets, and channel bit depths. A color lookup table (CLUT) 
component of the lutAtoBtype tag in ICC can be constructed as a 2x2x2 table derived from the 
matrix/offset/ID Luts CIELAB formula by running the eight corner input values (e.g., (0,0,0), 
(0,0,1), ... (1,1,1)) through the matrix. The non-canonical CIELAB encoding ranges and 
offsets can be rolled into the matrix/offset/ID Luts form, and the CLUT generated can thus 
handle the mapping of the non-canonical CIELAB encoding ranges and offsets. 
[0040] The CLUT entries can be made positive by subtracting the minimum CLUT entry 
from all and normalizing by dividing by the largest CLUT entry. The M curves can be three 
type-4 parametricCurveType tags of the form: 

(2) y = {ax + b) r + e for x>d and y = (cx + f) for x<d 

These curves can denormalize the output from the CLUT and apply the nonlinear function of 
the CIELAB to XYZ conversion, and scale by 32768/65535. 

[0041] A chromatic adaptation tag (chad) can be constructed for the illuminant specified in 
the non-canonical CIELAB parameterization, and this chad matrix can be used in the 
construction of the matrix stage of the lutAtoBtype tag in ICC. For example, a received white 
point can be expressed in degrees Kelvin, and this can be translated into an XYZ value for the 
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white point, which is used to calculate the chromatic adaptation matrix. The constructed matrix 
stage can denormalize the M-curve output and apply the adaptation from the parameterized 
CIELAB illuminant to the PCS D50 illuminant. Moreover, in the case of using the perceptual 
intent (one of the rendering intents available during gamut mapping), the constructed matrix 
stage for an AToB tag can map the parameterized CIELAB 's black point to the ICC perceptual 
intent black. 

[0042] The following description provides a derivation of JPX parameterized LAB profile 
elements to further illustrate the description above. Conversion of CIELAB to XYZ can be 
performed via the following equations: 



X = X w *f 
(3) Y = Y w *f 
Z = Z lv *f 



{L + 16) a 
116 + 500 



(Z + 16) " 
116 

"(1 + 16) b_ 

116 200 



o 

where X w , Y w , Z w are the tristimulous values of the white point, and f[x] = jc 3 for x > — , 

29 

m & f\A = "~ I x — — 1 for x < — . Equations (3) work for actual CIE L* a* b* units, but in 
841 ^ 29 J 29 

the case of JPX parameterized Lab data, the Lab values are encoded with specifiable ranges and 
offsets. Encoded JPX Lab can be decoded via the following equations: 

(4) a = — ^— (a -O a ) 
v/ 2" a — 1 



b = - 



2"> -1 



(K-o b ) 



where L e , a e , b e are the JPX encoded Lab values, R L , R a , R b are the JPX ranges for Lab (in 
Lab units), 0 L , O a , O b are the JPX Lab channel offsets, and n L , n a , n b are the JPX Lab 
original encoding bit depths. 
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[0043] In the context of an ICC v4 lutAtoB tag, the incoming L, a, b values ( L x , a i , b t ) are 

on [0..1] because they will already have been divided by (2 W - l), and thus, L i = — — — , 

2" L —1 



a e h 

a t = e - — b _ K 



L = R,L 



2 H -1 



, which allows equations (4) to be rewritten as: 



(5) a = R a a, a —O a 

b = R b b i —O b 

» > 2"" -1 4 

As the original JPX encoding ranges and bit depths are known, the offsets can be pre-scaled 
(0 L =-^ L T 1 °l> °a =^7Zi°<>> °» =^y^), (note that the offsets 0 L , O a , O b have 
now been redefined), and equations (5) can be simplified to: 

L = R L L i -0 L 

(6) a=R aai -O a 
b = R b b i -O b 



Modifying equations (3) above to reflect the conditions within the lutAtoB tag results in: 
X = X w *f 



(R L L i -0 L +\() + R a a i -O a 



(7) Y = Y w *f 
Z = Z w *f 



116 500 
(^Z,-Q, + 16) - 
116 

'(R^-O, +16) R b b t -O b 



116 



200 



[0044] Equations (3) above can be implemented as a matrix + offset followed by ID tables. 
Equations (3) can be rewritten as: 
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x=x w *f 

(8) Y = Y w *f 



116 



L + 



116 



500 



16 
116 



1 r 16 
L + 

116 116 



1 


( 


116 b 






L- 




V- 


116 


K 


200 j 


1 116 



The three terms that are each multiplied by in equations (8) above can be implemented as 
a 3x3 matrix: 



(9) 



B 



'U6\ 

0 
0 



where the 3x3matrix reduces to 



v 

1 

1 
1 



0 
0 

200 J) 



a 



' 116 ^ 

L + a 

500 
L 



200 



J 



125. 
0 



0 
0 



0 [-» 
I 50 



[0045] The rest of equations (8) can be captured in three ID Luts: 



(10) 



X w *f 



16 + A 

116 
16 + 5 



116 
16 + C 



116 



where A , B , C are the output of the matrix operation, and / is the same function used in 
equations (3). 

[0046] Equations (7) can be implemented in a similar manner with the addition of a vector 
of offsets to the form above: matrix + offsets vector, ID Lut form. Equations (7) can be 
rewritten as: 
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X = X w *f 
(11) Y = Y w *f 
2 = Z w *f 



116 



116 



ii6a 



R.L. + R a a i - 

' 500 " ' 500 J 



(16 -0 L ) 



116 



116 L ' 116 



1161, L ' 200 



200 



116 



FIG. 8 shows equations (1 1) with various terms highlighted. The boxed terms 800 can be 
implemented as a 3x3 matrix as follows: 



( r 



(12) 



B 



R, 



R r 



V 125 \ 
0 

0 



0 
0 

I 50 \ 



1 ' 500 

R L L i 

R^-^-R^ 
. 1 1 200 " 'J 



The circled terms 810 can be added as an offset vector: 

116 











(13) 


E 




B 




F 







( "60.] 




500 




0 




1160, 




[ 500 J 





R L L t + 



500 

R^-^RJ, 



R„a, - 



116*0 

500 
0 

1160, 



200 * ' 500 
The remaining circled terms 820 can be captured as three ID Luts: 



(14) 



f v\ 



x w 


*/ 


"16 


-0 L +D' 






116 






•/ 


"16 


-0 L +E~ 
116 




Z w 

V 


•/ 


"l6 


-0 L +F- 




H6 


) 



Note that 0 L is usually zero, as L* is positive. 

[0047] As the JPX parameterized Lab spaces can have varying white points, the conversion 
should go Lab->XYZjp X wp-»XYZp CS wp. The matrix of the lutAtoB tag can be used to handle 
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the XYZ white point conversion, the matrix + offset described above can be implemented in a 
2x2x2 CLUT in the lutAtoB tag, and the M-curves of the lutAtoB tag can handle the rest. 
[0048] A 2x2x2 CLUT that produces results equivalent to a 3x3 matrix can be generated by 
running the input corners of the CLUT through the matrix + offset. Using the matrix and offset 
from equations (12) and (13) above results in: 
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Rl 



V 



(15) 



0 
0 



Rl 


f 29 /? 
^125 


R, 

L 


0 


Rl 


0 


Rl 


f 29 ff 


Rl 


0 


«L 


0 






Rl 

L 


0 


Rl 


0 


Rl 




Rr 


0 




0 




<2L K 

^125 


"I 


0 


D 


A 

u 








0 


* £ 


0 


*i 


025 ° 


* £ 


0 




0 



0 
0 

^R b 

50 4 

0 
0 

29 „ ^ 











• 


0 












h 









+ 



29 V 



I, 125 
) 

O, 



0 

'29^ 



J 



50 



V 



V 



50 



J 



(f 



)) 

\ 



0 



R b 



o 

29 
'50 

0 
0 

50 * 

0 
0 

-R b 

50 * 



(0\ 
0 

ro^ 

i 



^ (f 



I 125 ") 
0 

(( 9Q 
- — 0 
I 125 "J 

e 0 ^ 



—R. o 



o 

50 



1 



rn 

0 



j 

29 Y\ 

a 1 



I, 125 " 
0 



V 

(f 



0 



50 \ 

l 125 ") 
0 

29 * 



V 

>7 



) 



0 

0 
29 

0 
0 

-^R>) 

50 ")) 



rn 
i 



v 1 ; 



v 



V 



125 
0 

1 50 4 J , 
I 125 'J 



125 
0 

'29 



\ 



O h 



J 



50 

' ' 29 o 

125 a 



V 



0 



(R a -o a ) 



125 



0 

50 * 



0 

19 (o 4 -«,)' 



V* 



V 50 

, 1 125 " 
' 29 

R L +T-o b 



50 



' 29 ^ 

R L -—o a 



125 



^ J) 

R L ^{R a -OS\ 

J 



125 



Rl 



29 v 



125 ° 



f 29 ^ 

['''so* J 

Rl 

R L +-(O b -R„)\ 
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This results in the following raw CLUT: 

< 125 j 



(16) 



000 
001 
010 

on 

100 
101 
110 

111 



125 
29 O a 
125 



'29 ^ 
' 29 

( ' 29 0 ^ 



R L - 



125 
29Q fl 
125 



29 S 

R L + — (R a -o a , 

L 125 V a a) I 



(o) 
(o) 
(o) 
(o) 



(29 0 b ^ 



29 
50 



50 



(O b -R b ) 



29 0^ 



50 



29 Oj^ 

\ ^ J 



(R L ) [Rl+^(0„ _R b ) 



( n 29 O b ^ 



V 



29 



(R L ) \R L +-(O b -R b ) 



[0049] However, an ICC v4 CLUT should have only positive values and produce output on 
[0.0.. 1 .0]. For a given raw CLUT, these conditions can be satisfied by determining minimum 
and maximum values, CLUTmin and CLUTmax, from the entries of the raw CLUT. All the 
values in the CLUT can then be forced to > 0 by subtracting CLUTmin from each entry, and 
output from the CLUT can be made to fall on [0.0.. 1.0] by dividing the elements of the raw 
CLUT by (CLUTmax - CLUTmin). This yields the following positive, normalized CLUT: 
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(17) 



290, 

— --CLUJmin 
125 



CLUTmz^CLUTmm 



290. 



125 



^--CLUTmin 



CLUTmax- CLUTmin 
f 29 



125 



(R a -O a -CLUTmin 



CLUTmax- CLUTmin 



f 29 
125 



(R a -O a -CLUTmin 



CLUTmaxr CLUTmin 
290 a 

R, --CLUTmin 

125 



CLUTmzx- CLUTmin 
290 



R, 



125 



--CLUTmin 



CLUTmax- CLUTmin 



R L +^ 5 (R a -O a )-CLUTmin 



CLUTmax- CLUTmin 



v 



J 



' R L +^R a -O a )-CLUTMn 
CL UTmax- CLUTmin 



VV 



-CLUTmin 



CLUTmax- CLUTmin 



-CLUTmin 
CLUTmax- CLUTmin. 



R L -CLUTmin 



CLUTmax- CLUTmin) 



R L - CLUTmin \ 
CLUTmax- CLUTmin) 



R L -CLUTmin 
CLUTmax- CLUTmin) 



' R L -CLUTmin > 
^CLUTmax- CLUTminj 



( 



\CL UTmax- CL C/2min, 



-CLUTmin "j 
CLUTmax- CLUTmin) 



-CLUTmin 



290 b 
50 



\ 



— CLUTmin 



CLUTmax- CLUTmin 

V 
'29 



50 



(O b -R„)- CLUTmin 



CLUTmax CLUTmin 



290, 



50 



-CLUTmin 



CLU'Enax- CLUTmin 



29 ^ 
— (O fc -&)- CZOTrnin 
50 



CLUTmax- CL UTmin 



290. 



/?, + ^-CIC/Knin 

* 50 



CLUTmax CLUTmin 



( 29 ^ 
R L +-(O b -R b )- CLUTmin 



CLUTmax- CLUTmin 



290. 



50 



-CLUTmin 



CZt/2rnax-CZ,t/7inin 



r 29 > 
R L +-{O b -R b )-CLUTmin 



CLUTmaxr CLUTmin 



J) 



Plugging the ICC PCS Lab Range and Offsets [R L = 100, 0 L = 0, Ra = 255, 0 a = 128, R b = 255, 
0b= 128] into (16) above results in the following raw CLUT: 
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(18) 



125 J 
3712S 

125 J 
f 3683^ 



V 125y 

f 3683^ 



V 125 y 
f 8788^ 



V 125y 

( 8788 N 



125 
^ 16183^1 



J 



125 ; 

16183 S 
125 J 



(0) 
(0) 
(0) 
(0) 
(100) 
(100) 
(100) 
(100) 



fl856> ' 

I 25 j 
3683 



v 50 
^ 1856' 



. 25 J 
' 3683^ 

» 50 J 
^ 4356 "| 

I 25 J 
(13171 
I 50 J 

( 1317 ^ 
{ 50 



'(-29.696) 
(-29.696) 
(29.464) 
(29.464) 
(70.304) 
(70.304) 
(129.464) 

Jl 29.464) 



(0.) 
(0.) 
(0.) 

(0.) 
(100.) 
(100.) 
(100.) 

(100.) 



(74.24) ^ 

(-73.66) 
(74.24) 
(-73.66) 
(174.24) 
(26.34) 
(174.24) 
(26.34) J 



which then gives a CLUTmin and CLUTmax as: CLUT min = - 



3683 
50 



= -73.66; 



CLUTmax = = 174.24 . 
25 

[0050J Subtracting by CLUTmin and dividing by (CLUTmax - CLUTmin) gives a positive 
normalized CLUT: 
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ft 



(19) 



10991 

<61975, 
( 10991' 

Ul975, 
|" 25781 N 

( 25781 s 

Ul975, 
35991 N 

.61975, 
f 35991 N 

Ul975, 
5078 1 N 

.61975, 

( 5on\\ 

J 



f 



3683 
^12395 
f 3683 

^12395 
( 3683 

U2395 y 
( 3683 ^ 

V12395; 
( 8683 \ 
U2395 
( 8683 
U2395 y 
( 8683 \ 
U2395 



1479 



7 



8683 ^ 



^61975 J U2395J U479 



U479J 
(0) 

( 1479 ^1 
U479J 

(0) 

(1) 
( 1000 

U479 

(1) 

( 1000^ 



J) 



'(0.177346) 
(0.177346) 
(0.41599) 
(0.41599) 
(0.580734) 
(0.580734) 
(0.819379) 
(0.819379) 



(0.291736) 
(0.297136) 
(0.297136) 
(0.297136) 
(0.700524) 
(0.700524) 
(0.700524) 
(0.700524) 



(0.596612)") 

(0.) 
(0.596612) 

(0.) 

(1.) 
(0.403388) 

(1) 
(0.403388) 



Multiplying this by 65535 gives the following hex entries for the ICC PCS Lab CLUT: 



(20) 



'(2<*66 16 ) 


(4cll 16 ) 


(9Sbb l6 )) 


(2</66 16 ) 


(4cll 16 ) 


(0,.) 


(6a7e l6 ) 


(4cll 16 ) 


(9«Mi«) 


(6a7e l6 ) 


(4cll 16 ) 


(0i«) 


(94aa 16 ) 


(M55 16 ) 




(94aa 16 ) 


(M55 16 ) 


(6744 16 ) 


(dlc2 l6 ) 


(M55 16 ) 


(Mis) 


^2 16 ) 


(M55 16 ) 


(6744 16 )J 



This CLUT thus handles the processing up through equation (13) above, and the M curves can 
handle equation (14). 

[0051] As a check regarding adding the offset to each of the CLUT corners, as in equations 
(15), consider bilinear interpolation: the offset after interpolation, x = (ax { + bx 2 ) + O x , is not 
equal to the offset before interpolation, x = a{x x + O x ) + b (x 2 + O x ) ; but for linear 
interpolation with b = (1 — a) , simplifying these equations results in: 
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(21) 



x = [a(x t + (1 - a)x 2 ) + O x ] = ax l +O x -(-l + a)x 2 

x = [a(x { + O x ) + (1 - a)(x 2 + O x )] = O x + ax x - (-1 + a)x 2 



Thus, they are equivalent. Here is a full comparison of matrix + offsets with CLUT with offsets 
added before interpolation: 



(22) 



e \ e 2 e 3 



e 4 e S e 6 



'10 



\ e nj 



Le, + ae-, + be, + e, n ^ 



Le 4 + ae 5 + be 6 + e, 
Le 7 + ae % + be 6 + e ]2 j 



Doing this via CLUT with offsets added to the CLUT: 



(23) ml = 



1 ' e \ e 2 

e 4 e 5 e 6 



Matrix Form [CLUT1 = 



{ml. {{0}, 

ml. {{0}, 

ml. {{0}, 

ml. {{0}, 

ml. {{1}, 

ml. {{1}, 

ml. {{I}, 

ml. {{1}, 



{0}, {0}}, 

{0}, {1}}, 

{1}, {0}}, 

0}, «}, 

{0}, {0}}, 

{0} {1}}, 

{1}, {0}}, 

{1}. {1}} 



}] 
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(0) 


(0) 


(0) 




M 

v o / 


\ y J 




(eA 




(e 2 +e 3 ) 


(e 5 +e 6 ) 


(e s +e 9 ) 




(e 4 ) 


(e 7 ) 


(e ] +e 3 ) 


(e A +e 6 ) 


(e 7 + <? 9 ) 


(e, +e 2 ) 


(e 4 +e 5 ) 


(e 7 +e 8 ) 



{(e { +e 2 +e 3 ) (e 4 + e 5 +e 6 ) (e 7 + e 8 + e 9 )^ 

Adding in the offsets: 

(24) MatrixFormfCLUTwOffs = Map[# + {elO, el 1, el2} &, CLUT1]] 



(ej 

( e 2 +e i0 ) 

(e 2 +e 3 +e }0 ) 

( e !+ e io) 
(e,+e 3 +e 10 ) 

(e,+e 2 +e 10 ) 



( e 6+0 
(e 5 +<? u ) 
(e 5 +e 6 +e u ) 

(e 4 +e u ) 
(e 4 +e 6 +e u ) 
(e 4 +e 5 +e n ) 



CO 
(e 9 +e n ) 

(e s +e 9 + e l2 ) 

(e 7 +e u ) 
(e 7 +e 9 +e i2 ) 
(e 7 +e 8 +e 12 ) 



(e, +e 2 +e 3 +e l0 ) (e 4 +e 5 +e 6 +e u ) (e 7 +e 8 +e 9 +e l2 \ 



And defining a trilinear interpolator: 
trilin[CLUT_, Xl_, X2_, X3J:= 

Module [{AO, Al, A2, A3, A4, A5, A6, A7}, 
AO - CLUT [[1]]; 
Al = CLUT [[5]] - CLUT [[1]]; 

A2 = CLUT [[3]] - CLUT [[1]]; 

A3 = CLUT [[2]] - CLUT [[1]]; 

A4 = CLUT [[7]] - CLUT [[5]] - CLUT [[3]] + CLUT [[1]]; 
A5 = CLUT [[6]] - CLUT [[2]] - CLUT [[5]] + CLUT [[1]]; 
A6 = CLUT [[4]] - CLUT [[2]] - CLUT [[3]] + CLUT [[1]]; 
A7 = CLUT [[8]] - CLUT [[4]] - CLUT [[7]] + CLUT [[5]] + 
CLUT [[2]] + 
CLUT [[3]] - CLUT [[1]] - CLUT [[6]]; 
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(AO + (Al * XI) + (A2 * X2) + (A3 * X3) + (XI * A2) + 

A5 (XI * X3) + A6 (X2 * X3) + A7 (XI * X2 / X3)) 

]; 



MatrixForm[trilin[CLUTwOffs, L, a, b]] 



t Le x + ae 2 + be 3 + e l0 ^ 



Le 4 +ae 5 be 6 +e u 
^Le 7 + ae g + be 9 + e n j 

which equals (22) above. 

[0052] With respect to the M curves, equation (14) is now presented again with the 0 L 
dropped, as 0 L is zero in this implementation, and D, E, F are the output from the CLUT 









x w 


*f 


(25) 


Y 






*f 








z w 


*f 



16 + D 

116 
16 + E 



116 

16 + F 



116 



The output of the CLUT is on [0.0.. 1.0], and the denormalized output value can be obtained by 
multiplying by (CLUTmax-CLUTmin) and then adding CLUTmin: 



(26) 



f Y\ 



x w *f 

Y w *f 
Zw*f 



16 + DjCLUTmzK- CLUTmin) + CLUTmin 
16 + E(CLUT max- CLUTmin) + CLUTmin 



116 

16 + F(CLUT max- CLUTmin) + CLUTmin 



116 



Looking at the function / again, f[x] = x 3 for x > ^ , and f[x] = 



4 . 6 

x for x < — , 

29 J 29 



and plugging in the terms for D from equation (26) into the function / results in: 
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/ 

(27) 



16 + D(CLUT max- CLUT min) 



+ CLUT min ~j_ 



p6 + 



116 

.3 



D{CLUT max- CLUTmin) + CLUT min 
116 



. 16 + D(CLUT max- CLUT min) + CLUT min 6 

lor > — , which resolves to 

116 29 

CLUT min- 8 

£) > , and 

CL UT max- CLUT min 



/ 



16 + £>(CXE/r max- CZLTmin) + CELT min 
116 



(28) 108 { 16 + £>(C£LTmax- CELT min)+ CZLTmin 4 ^| 
841 { " 116 29 J ~ 

27(Z)(CZLTmax- CLUTmin)+ CLUTmin) 
24389 

, \6 + D{CLUT max- CLUT min) -CLUT min 6 .. . 

lor < — , which resolves to 

116 29 

CZLTmin-8 

D< — 



CLUT max- CLUT 'min 
[0053] This is done as follows: 

0 , r 16 + D(CZLTmax-Cmmin+CI,LTmin 6 „ 

Solve [ = — ,D] 

116 29 

ff _ -8 + CLUTmin 

{{£>-> }} 

CZLTmax- CZLTmin 

32768 

[0054] Then by scaling the output from the M curves by to get it onto the correct 

65535 

for the matrix stage, the equations in (26) become: 

(29) X = X { 32768 Y 16 + -P(CELTmax- CLUT min) + CZLTmin V 
"\ 65535 A 116 

. _^ CZLTmin-8 
for D > 

CL UT max- CL UT min 

and 
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X = X, 



32768^1 21(D(CLUT max- CLUT min) + CLUT min) 



^65535 J 
for 73 < - 



24389 



CLUT min- 8 



CLUTmax- CLUTmin 



Y = y f 32768 Y 16 + E(CLUT max- CLUTmin) + CLUTmin 
"\ 65535 A 116 

for E> — 



CLUT min- 8 



CZC/r max- CLUT min 



and 



r = r. 



w 



32768 



U5535J. 
for £ < - 



21(E(CLU1 'max- CIE/r min) + CXE/r min) 
24389 

CLUT min- 8 



Z = 2 



CLUTmax- CLUTmin 

32768 Y 16 + F(CLUT max- CLUTmi n) + CLUTmin ^ 

J 



.65535 J 1 
for F > - 



24389 



CLUT min- 8 



CZ£/r max- C£t/r min 



and 



Z = Z 



32768^1 21(F(CLUT max- CZ f/J min) + CL UT min 



65535 J" 
for F < — 



24389 



CLLT min- 8 



CX£/r max- CLt/7/ min 
[0055] The equations (29) can be represented as ICC type 4 parametric curves, which have the 
following structure: y = (ax + b) r + e for (x > d) , and y = (cx + /) for (x<d). Looking at just 
the X cases in (29), the initial factor can be moved inside the cubed expression by taking its 
cube root, and the cubed portion in (29) can be rewritten to result in: 



(30) X = 



X 



( 32768 Y) * ( (CLUTmax- CLUTmin) n 16 + CLUTmin 



%65535jJ I 



116 



116 
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for D> — 



CLUT min-8 



CLUT max- CLUT min 
which puts this in a form that matches y = (ax + b) r + e for x > d : 

X = 



(31) 



32768 Y) 3 ( (CLUT max- CLUT min) N 



65535 J J { 



116 



( v ( 32768 \)~ 3 ( 16 + CLUT min)} 3 , _^ CLUT min- 8 

X w\ 77777 I 777 I / or D * ~ 



1,65535 



116 



CLUT max- CLUT min 



resulting in: 



y = 3 



a = 



b = 



( (CLUT max- CLUT min) Y ( 32768 



116 



16 + C£ £77^^(732768 



65535 J 



d = — 



116 
CZt/r min- 8 



J, 65535. 



■ w 



CLUT max- CLUT min 



e = 0 

[0056] For the second form for X in (29): 



X = X ( 32768 1 W (D(CLUT max- CLUT min) + CLUT min) 
*t 65535 J 24389 



for D< — 



CZ IT min- 8 



CZ C/T max- CL UT min 
which can be put into y = cx + f form: 

(32) x -- x -W^^ CLUT ™- cLuTm[n)D+ 



(" 32768 ^ 27C/XT7' min 
9 1,65535 J 24389 
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. _ CLUT min-S 

for D<-- 



CLUT max- CLUT min 
which simplifies to: 

294912 

X = (CLUT max- CLUT min)X lv D + 

532777705 w 

294912 _____ . . _ CLUT mm-Z 

CLUT mm X w for £)<-- 



532777705 CZf/r max- min 



and this results in: 



294912 

c = (CLUT max- CLUT min)__ „, 

532777705 * 

294912 

/ = CLUT min JT„ 

632777705 r 

[0057] Here are Mathematica® simplifications of the two terms in (32): 

( 32768Y 27 ^ 

Simplify [X„ (CLUTmax - CLUTmin)D] 

\65535 )\ 24389 J 

2949\2(CLUT max- CLUT min)DX w 



532777705 



Simplify [X w 



( 32768 



1,65535 
294912CZf/rminX 



27 'CLUT mm 



24389 



] 



532777705 

So the parameters for the type 4 parameteric curve tags are : 
y = 3 



294912 , , 

c = (CLUT max- CL UT mm)[X w \Y W \Z w ] 

532777705 r| r| w 

CLUT min- 8 

a = -- 



CZt/T max- CL_/T min 
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e = 0 



294912 i , 

/ = CLUT mm[X w \Y W \Z w ] 

532777705 r| w 

[0058] A routine to evaluate these parameters for given CLUTmax, CLUTmin, and white point 

XYZ can be as follows: 

McurveFunc[CLUTmax_, CLUTmin_, {X_, Y_, Z_}]:= 
Module[{a, b, c, d, e, f}, 



a = { 



( (CLUT max- CLUT min) 



V 



116 



732768 
1,65535. 



X 



(CLUT max = CLUT min) Y ( 32768 V 
116 XI 65535 J 

(CL UT max- CL UT min) Y ( 32768 ^| 
116 XU5535/ 



Print["a=", N[a]]; 



b = 



16 + CLUT min Yf 32768^1 



65535 



116 

r l_6 + CLt7JminY f 32768^1 

v 



\ 



X 

j ) 



116 



65535 



AV 

'16 + CZt/rminYf 32768 



v 



116 



J 



V 



65535 



Print["b=',N[b]]; 
294912 

c = {77^7777 (CLC/r max- CLUT min)(X), 



532777705 
294912 



532777705 
294912 



(CZC/7 max- CZt/r min)(y), 
(CIC/r max- CZ £77/ min)(Z)} ; 



532777705 
Print["c=", N[c]]; 
d _ CLUT min-S 

~ CLUT max- CLUT min 
Print["d=", N[d]]; 
Print["e=", O]; 
29491 2 

/ = { — — CLUT mm(x\ 



29 



Patent Application 

Att. Dkt. 07844-647001 / P600 



294912 

532777705 
294912 



CLUTmm(Y\ 
CLUTmin(z)}; 



] 



532777705 
Print["£=", N[f}}; 
Print[N[{a,b,c,d,{0,0,0},f}]] ; 
Print["in hex:"]; 
Print[ 

BaseForm[MatrixForm[Round[{a, b, c, d,{0,0,0},f * 65536.0]], 16]]; 
Print["dec from hex:"]; 
Print[ 

MatrixForm[N[Round[{a, b, c, d{0,0,0}f} * 65536.0]/65536.0]]]; 



[0059] For the ICC PCS D50 Lab case: 

McurveFunc [174.24, -73.66, {0.9642, 1.0, 0.8249}] 

a= {1.67571, 1.6962, 1.59078} 

b = {-0.389761, -0.394526, -0.370006} 

c= {0.132309, 0.137222, 0.113194} 

d= {0.329407} 

e = 0 

f = [-0.0393138, -0.0407735, -0.0336341} 

{{1.67571, 1.6962, 1.59078}, 
{-0.389761, -0.394526, -0.370006}, 
{0.132309, 0.137222, 0.113194}, {0.329407}, 
{0., 0., 0.}, {-0.0393138, -0.0407735, -0.0336341}} 

in hex: 

' {/ac/c I6 ,l&23a 16 ,1973e I6 } > 

{-63c7 16 ,-6500 16 ,-5e69 16 } 

{2\df x 2321, lcfaj 
{5454 I6 } 

{0,6,0, 6 A 6 } 

{- fl 10 16 ,-fl70 16 ,-89c l6 } 
dec from hex: 

{1.67572,16962,1.59079} 
{- 0.389755,-0.39453 1,-0.37001} 
{0.132309,0.137222,0.11319} 
{0.329407} 
{0.,0.,0.} 

J- 0.0393066 -0.0407715 -0.0336304}, 
negative hex terms: 
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- 63c7 16 = 0xFFFF9C39 ,-6500 16 = 0xFFFF9B00, 

- 5eb9 x6 = 0xFFFFF5F0, 
-al0 X6 =0xFFFFF5F0, 

-a70 16 = OxFFFFF 590,-S9c l6 = 0xFFFFF764 

where, y = 3.0 = 0x00030000 . 

[0060] For the A2B1 tag, the matrix is an identity matrix with 0,0,0 offsets. For the A2B0 tag, 
the matrix is non-identity as it handles the perceptual black point. If this profile has a D50 
white point, there is no chromatic adaptation to be handled in the matrices. Here's the 
formulation for the matrix to handle the perceptual black point: 



((( 



e 2 e 3 



S S e 9 J 



""10 



\ e nJ 



-» 



• bperc 



L Wpcs j 



0 



0 



0 

( Y > 

Y bperc 

. ~~Y , 

y wpcs J 



vv 



1— 



0 
0 

z 



bperc 



'wpcs J J 



.chad 



L bperc 



bperc 



with 



fx > 

wpcs 




( 0.9642^ 


Y 

wpcs 




1.0000 


z 

\ wpcs J 




^0.8249^ 



X, 



bperc 



and 

" ( 0.00336 > 



Z 

\ bperc J 



0.0034731 
0.00287 



J 



giving 
^0.996515 



0 
0 



0 

0.996527 
0 



0 
0 

0.996521 







( 0.00336 > 


.chad 


+ 


0.0034731 


) J 




^ 0.00287 j 



As this profile for the ICC PCS has a D50 white point, the chad is an identify matrix: 



'0.996515 


0 


0 > 






0 


0^ 




( 0.00336 > 


0 


0.996527 


0 


• 


0 


1 


0 


• 


0.0034731 


, o 


0 


.996521; 




,o 


0 


h 




^ 0.00287 ; 



so the matrix and offsets are: 
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'0.9965 15 0 
0 0.996527 



v 



0 W 0.00336 ^ 
0 + 0.0034731 
0 0 0.996521 J { 0.00287 

and in hex this is: 

'0000FF15 00000000 00000000^ 

00000000 0000FF1C 00000000 
^ 00000000 00000000 OOOOFFW; 
[0061] For the A2B1 tag, there is no perceptual black point to account to deal with, so the 
matrix and offsets can be: 



fOOOOOODC^ 
000000F1 
OOOOOOtfC 



( 1 


0 


°1 






0 


1 


0 


+ 


0 


,0 


0 


0 




,0; 



in hex: 



00010000 00000000 00000000 
00000000 00010000 00000000 
00000000 00000000 00010000 



+ 



00000000^ 

00000000 

00000000 



A-Curves and B-Curves can both be identity curveType tags, identified by a 0 count field. 
[0062] FIG. 9 illustrates creation of a transform-defining element according to another 
implementation. A first processing stage that accounts for the range of a parameterized color 
space is determined at 900. A second processing stage that accounts for the offset of the 
parameterized color space and defines at least a portion of a conversion of the image color 
space to a chromatic adaptation color space is determined at 910. A third processing stage that 
defines a chromatic adaptation according to the white point of the parameterized color space is 
determined at 920. The transform-defining element illustrated in FIG. 7 can be used. Thus, 
mapping of non-canonical ranges and channel bit depths can be defined in a three dimensional 
interpolation table stage of the transform-defining element 450, mapping of non-canonical 
offsets can be defined in the one dimensional transforms, and non-canonical white points can be 
compensated for in a matrix stage of the transform-defining element 450. 
[0063] In the context of constructing an ICC v4 profile with a CIEXYZ PCS, the invention can 
combine the conversion from a non-canonical CEELAB encoding with the conversion from 
CIELAB to the CIEXYZ D50 PCS, as follows, resulting in color profiles that are independent 
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of the data channel bit depths. In an ICC v4 profile with a CIEXYZ PCS, the lutAtoBtags can 
be constructed to implement the CIELAB to XYZ conversion as above, but with non-canonical 
offsets being delayed into the M curves rather than being applied in the matrix used in 
constructing the CLUT. The input and output of the CLUT can use the full range, 0..1, for 
maximum precision. This results in improved accuracy. 

[0064] The CLUT entries can be made positive by negating the incoming CIELAB B channel 
(1-B) and swapping the rows in the CLUT that have a 1 in the B channel with those that have a 
0 in the B channel. As before, the M curves can be three type-4 parametricCurveType tags of 
the form: 

(33) y = (ax + b) r +e for x>d and y = (cx + /) for x<d 

These curves can denormalize the output from the CLUT, apply the nonlinear function of the 
CIELAB to XYZ conversion, and apply the non-canonical channel offsets. The input and 
output of the M curves can use the full range, 0..1, for maximum precision. By arranging the 
operations appropriately, the data used in making the CLUT can be more linearly related, and 
thus linear interpolation with the CLUT produces minimal interpolation errors. 
[0065] Additionally, a chad can be constructed for the illuminant specified in the non-canonical 
CIELAB parameterization, and this chad matrix can be used in the construction of the matrix 
stage of the lutAtoBtype tag in ICC. The constructed matrix stage can denormalize the M- 
curve output, multiply by the LAB space white point XYZ, scale by 32768/65535, and apply 
the adaptation from the parameterized CIELAB illuminant to the PCS D50 illuminant. 
Moreover, in the case of using the perceptual intent (one of the rendering intents available 
during gamut mappint), the constructed matrix stage for an AToB tag can map the 
parameterized CIELAB ? s black point to the ICC perceptual intent black. 
[0066] The following description provides a derivation of JPX parameterized LAB profile 
elements to further illustrate the description above. Conversion of CIELAB to XYZ can be 
performed via the following equations: 



X = X w *f 



(£ + 16) a 
116 + 500 
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Y = Y w *f 
Z = Z w *f 



{L + 16) 

116 
1 + 16 



116 200 



where X W ,Y W , Z w are the tristimulous values of the white point and f[x] = x* for x > — , 

29 



108 



841 



x — 



29 
ForX 



for x < — . In term of the current ranges and offsets, these become: 
29 



X = X w *f 



Rl (A» - O. „ ) + 1 6 R a (a.-O n 

L \ in Ln / _|_ a v m <? 



116 500 
29 116 V n/ 500 V 



where: 

A* » a « > *e jpx encoded Lab values 

R L ,R a , R b are the jpx ranges for Lab (note that these are in Lab units) 
O l ,O q , O b are the jpx Lab channel offsets 
n L > n a > n b ^ JP X Lab or original encoding bit - depths 
normalized offsets 

o = ° L 
O 

o - . °» . 

*• ( 2 "*-l) 



[0067] This can be rearranged to: 
X = X w f 



R L r R a 

— + — s -a. + 
116 " 500 " 



29 116 



500 



in terms of the current notation: 



X = X w f 



L in a in 



{29 



-s L o Ln -s a o a 



where ^ S„ S 6 S, =5, +5 fl , 5, =5 i5 5 2 =5, + 

O x = O to * 5, + O an *S a ,O y = 0 Ln *S L ,0 2 = 0 Ln *S L +(\-O bn )*S b . 
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[0068] ForY: 



Y = Y w *f 



116 



which can be rearranged to: 



Y = Y„*f 



116 



{29 116 L \ 



in terms of the current notation: 



Y = Y w *f 
[0069] ForZ: 



5^+1—-^ 



116 



this can be rearranged to 

r R l 



Z = Z w *f 
further to: 

z = z„*/ 



200 
16 & 



-—o Ln +— -^-k +-*-a 



116 " 116 



116 200 



200 



6n 



116 " 200 " 
in terais of the current notation: 



- — + f — - — 0, + ^- O bn 
{29 116 Ln 200 



S L L in- S b b in + 



[0070] These 3 equations result in the following matrix and offset form: 
X = X w f 



Z = Z w *f 



(4 > 
S,L in +S a a. + SrO^-SO 

l in a in <^c\ L Ln a an 

v2y )\ 

S L L in +[^-S L 0 Ln 
S L L in -S b b in+ (^--S L 0 Ln+ S b O b ^ 



{29 



0 



matrix: 

(S L S a 

s L o 0 

A 0 

If any b input is replaced with (1-b), the result is a positive matrix: 
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S L 0 

k ° 



0 



normalized matrix: 



i 0 



0 

5, 



v5 z 



Offsets vector: 
' 4 

-^+o btt s b -o Ln s L 

[0071] This matrix gives the following CLUT: 



in 



fo o 6) 

0 0 1 

0 1 0 

0 1 1 

1 0 0 
1 0 1 
1 1 0 

1 1 1 



CLUT 






0 


0 


0 

s„ 


0 


0 




s a 






0 


0 


s x 






s 0 


0 


s b 


s x 


s 2 


s L 


s L 


s L 


s x 


s, 




s t 


s L 


s L +s„ 


s x 




s, 


s L +s a 


s L 


s L 


s x 


s y 


s 2 


s L +s a 


s L 


s L +s b 


s x 







[0072] Using (1-b) instead of b can be compensated for by swapping rows with zero and one in 
the third input column: 
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(0 0 0) 

0 0 1 

0 1 C 

0 1 1 



0 0 
0 1 



0 



0 



s 

i 

s x 
s± 

s x 

s± 

s x 
s L + s a 

s x 
s L +s a 



0 

0 
0 



0 



0 



Sl_ s l +s„ 



s L + s b 
s 2 

Sz 



Noting that -f- = 1, S L +S a = S x , and S L + S b = S z , the resulting CLUT is: 



(0 0 0^ 
0 0 1 



1 0 0 



0 

1 



1 1 



with offsets: 

( A 



0 0-^ 



0 

s x 

s x 

s x 

s x 
1 

1 



0 

0 
0 

0 

1 

1 

1 
1 



0 

s z 

0 

1 

s z 
1 

s 



2 _ 



—-o an s a -o Ln s L 
—-o Ln s, 



29 



+ o bn s b -o Ln s L 



[0073] Recalling that 
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o x =o Ln *s L+ o an *s a 
o y =o Ln *s L 

o 2 =o Ul *s L+ (i-o bn )*s b 

-*o 2 =o Ln *s L -o bn s b +s b 

the offset vector can be rewritten as: 
f 4 



29 



I 29 



-(o an s a +o Ln s L ) 
-o Ln s L -o bn s b 



substituting the O xyz (note: {0 Ln S L - O bn S„ ) = 0 2 -S b ): 



29 * 

giving the offset vector as: 
f 4 



0\ = 
0' = 



-o x 

29 x 
T9-°> 

4 

— + S b -0 2 
29 " 2 



[0074] The M-curves can denormalize the output of the CLUT, apply the offsets 0' xyz , and 
apply the function f from the Lab to XYZ conversion: 

f[x] = x 3 for x> — 
J 29 



/[*] = 



108 
841 



4 V 6 
for x < — 
29 



29 



As before, the M-curves can be type 4 parametric curves of the form: y = (ax + b) r + e for 
(x > d) , and y = (cx + f) for (x < d) . The general form for the m-curves can be: 



If 



, 29 J 
U=T 3 



else 
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" 8411, " 29 J 

[0075] The m-curves can be normalized so that output is on [0.. 1] by dividing by the max value 
of U n : max U n = r„ 3 = (S n + 0\ ) 3 . So for the channel 1 curve: 



T x =X*S x +0' x 
k[(X*S x +0' x )>-^ 



->If 



x*s x > 



-a 



►if 



x> 



r 6_ 
1 29 

29 



x ~ (s x +o< x y 



X*S r 0\ 



K S X +0' X S x +0' x ; 



else 



U = 



108 /" 
841 1 



X*S x+ 0' x - 



29, 



108 _ 108 
) 841 J , 841 



* 29, 



' (s x+ o' x y (s x+ o> x y 

[0076] This results in type 4 parametric curve parameters as follows: 



a. = 



S„+0\ 



b = 



0\ 



c. = 



S x+ 0' x 



108 

84 r 



29 



— + S x -O x 
29 

~-o, 

29 x 

4 

—+s-o x 

29 x x 



108 
841 x 



29 + S " ° x . 



_6_ 
29 



fi-V 

V29 ' 



29 x 



e x =0 
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— {a -— 
(s x +o\) 



simialarily for channel 2: 

Yy =3 



108 
841 



U-o)-± 

{{29 x ) 29 



108 
841 



—+s x -o x 

29 



— +s -o, 

29 X X 
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— + Sy -Oy 

29 r Y 
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+ Sy~O y 
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108 
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J v 29 
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841 1 r 29, 
(5 y+ OV) 



and for channel 3: 
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_ 29 * 



_6_ 
29 



e z =0 



/z = 
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841 



rr 4 c 
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[0077] For the matrix stage: 





e 2 






^10 








+ 






e 8 


e 9> 







The matrix denormalizes the output from the m-curves, multiplies by XYZ W and scales by 
32768 

— — , and applies the chromatic adaptation via the chad: rel col matrix = chad.scalingMatrix, 
65535 

where the Scaling Matrix is: 



655535 \ 29 x x 
0 



0 



0 

65535 "\ 29 Y Y 



0 



0 
0 



65535 T29 



The relative colorimetric offset vector = 



fo^ 

0 



Perceptual Matrix = PercBPScalingMatrix.rel col matrix. 
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PercBPScalingMatrix = 



1-- 



bperc 



0 



' wpcs J 



0 



0 



1- 



bperc 



wpcs J 
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1— 



0 
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bperc 



wpcs J J 



Perceptual offset vector = 



32768 



' bperc 



65535 
32768 

65535 bperc 
32768 

65535 *" erc 



with 



wpcs 



wpcs 

z 

\ wpcs J 



( 0.9642^ 
1.0000 
0.8249, 
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fx } 

A bperc 




( 32768^ 




^ 0.00336 ' 




( 32768 \_ 


'0.00168003^ 
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bperc 
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0.0034731 
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0.00173658 




[ 65535 J 




^ 0.00287 J 
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bperc J 








v 0.00143502 y 



[0078] The invention and all of the functional operations described in this specification can be 
implemented in digital electronic circuitry, or in computer hardware, firmware, software, or in 
combinations of them, such as the structural means disclosed in this specification and structural 
equivalents thereof. Apparatus of the invention can be implemented in a software product (e.g., 
a computer program product) tangibly embodied in a machine-readable storage device for 
execution by a programmable processor; and processing operations of the invention can be 
performed by a programmable processor executing a program of instructions to perform 
functions of the invention by operating on input data and generating output. The invention can 
be implemented advantageously in one or more software programs that are executable on a 
programmable system including at least one programmable processor coupled to receive data 
and instructions from, and to transmit data and instructions to, a data storage system, at least 
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one input device, and at least one output device. Each software program can be implemented in 
a high-level procedural or object-oriented programming language, or in assembly or machine 
language if desired; and in any case, the language can be a compiled or interpreted language. 
Suitable processors include, by way of example, both general and special purpose 
microprocessors. Generally, a processor will receive instructions and data from a read-only 
memory, a random access memory and/or a machine-readable signal (e.g., a digital signal 
received through a network connection). Generally, a computer will include one or more mass 
storage devices for storing data files; such devices include magnetic disks, such as internal hard 
disks and removable disks, magneto-optical disks, and optical disks. Storage devices suitable 
for tangibly embodying software program instructions and data include all forms of non-volatile 
memory, including by way of example semiconductor memory devices, such as EPROM 
(electrically programmable read-only memory), EEPROM (electrically erasable programmable 
read-only memory), and flash memory devices; magnetic disks such as internal hard disks and 
removable disks; magneto-optical disks; and optical disks, such as CD-ROM disks;. Any of the 
foregoing can be supplemented by, or incorporated in, ASICs (application-specific integrated 
circuits). 

[0079] To provide for interaction with a user, the invention can be implemented on a computer 
system having a display device such as a monitor or LCD (liquid crystal display) screen for 
displaying information to the user and a keyboard and a pointing device such as a mouse or a 
trackball by which the user can provide input to the computer system. The computer system 
can be programmed to provide a graphical user interface through which computer programs 
interact with users. 

[0080] The invention has been described in terms of particular embodiments. Other 
embodiments are within the scope of the following claims. For example, the systems and 
techniques described can be applied to other color spaces, such as the YCC (Luma-Chroma- 
Chroma) color space. Moreover, other defined color profile architectures and transform- 
defining elements can be used, such as a CSA (color space array) and a CRD (color rendering 
dictionary) in the PostScript® language. The operations of the invention can be performed in a 
different order and still achieve desirable results. 
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